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Abstract 

A phenomenological theory of mode collapse-revival for a system with a weak 
nonlinearity is presented. The theory takes into account fluctuations of the 
number n of quanta of oscillations. The collapse time r^ 1 for a mode of 
frequency ujq turns out to be uJo(y/n/2hujo \ v |), where v is a coefficient in 
a nonlinear correction to the mode frequency with respect to the oscillation 
energy E. For a Bose gas in a harmonic trap r c is order of uJq 1 \J fiN/EhujQ. 
This value is order of 250 ms. for typical experimental conditions. The 
coefficient v is calculated for the breathing mode in an isotropic trap. 
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The recent discovery of Bose-Einstein condensation of alkali atoms |IJ confined in mag- 
netic traps has opened a new important field of investigation of quantum phenomena on 
a macroscopic scale. The most interesting possibility is to observe specific quantum phe- 
nomena which have no classical analogy. One such phenomenon is the collapse and revival 
of coherent quantum states recently observed for an atom in a electromagnetic field in a 
resonant cavity |§ . The problem has a developed theoretical background, and for a com- 
prehensive discussion see in Ref. ||. 

Several significant works have been carried out recently about manifestation of this phe- 
nomenon in Bose-condensed gases. In paper 0] properties of a ground state are discussed. 
Kuklov et. al. || considered dephasing of a finite amplitude oscillation in a finite Bose 
system because of atom occupation number fluctuations in the oscillation. They suggested 
that this mechanism can explain the observed damping of oscillations and used a microscopic 
model to calculate the collapse time and they concluded that the theory can give a reason- 
able order of magnitude of this time compared to the damping observed in experiments ||, 
0. In [|] the authors considered the dephasing because of the zero-point fluctuations of 
occupation numbers of atoms. Since the authors have withdrawn the paper we don't discuss 
it in detail. Note only that we believe that zero-point fluctuations cannot result in the effect 
under consideration. 

In this paper we will present a general phenomenological theory of collapse-revival of a 
highly excited mode of oscillation of a superfluid Bose system. We will take into account not 
fluctuations of occupation numbers of atoms but rather fluctuations in the numbers of quanta 
of oscillations of the mode. (This is by no means the same for low-frequency oscillations.) We 
believe that this choice of variables is more proper then used in ||. The physical meaning 
of the effect is that a coherent oscillation is a linear combination of stationary states of 
a system and the number of quanta is a natural specification of a stationary state of an 
oscillator. In this way one doesn't meet the troublesome problem of diagonalization of the 
Hamiltonian. Of course a problem can be solved in any variables and a microscopic theory 
always gives more information about some mechanism of the effect. From another side, the 
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phenomenological approach is more transparent and gives direct possibility for comparison 
with an experiment. 

According to general theory [Q] the effect in consideration is based on the fact that 
transition frequencies for stationary component of a coherent state of the system are different. 
This effect doesn't exist for a purely harmonic oscillator and hence is completely defined 
by non-linearity. In § a two-level system in an electromagnetic field has been considered. 
One may say that this is a limiting case of an extremely strong nonlinearity. For our system 
however the opposite case of weak nonlinearity is more proper. 

First of all we will present the collapse-revival theory in a form more suitable for our 
purpose. In the case of weak nonlinearity the frequency of the oscillation mode under 
consideration can be written as: 

UJ = UJq + 5u = u (l + vE), (1) 

where E is an energy of oscillations, and | v \ E « 1 (weak non-linearity). Stress that 
constant v describing the first nonlinear correction has a purely classical meaning. For an 
anharmonic oscillator with a Hamiltonian H = p 2 /2m + muj\x 2 : /2 + ax 3 /3 + j3x 4 /4 one has: 

5a 2 3/3 
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Taking into account that hu = (dE n /dn), we can rewrite this equation in the quantum form 

uj n = E n /h = u n + bn 2 /2, (3) 

where n is the number of quanta in the given excited state of the system (n » 1), E n is 
the energy of the state and 

b = hujQis. (4) 

Let us consider now a coherent state of our oscillator. Its wave function can be presented 
in the well-known form: 

ip = C^2c n i[) n exp(-iuj n ), | c n \ 2 = exp(-n)n n /n\, (5) 



where ip n is a wave function of a stationary state of the oscillator where n quanta are exited, 
fiuiQfi = E. All quantities besides u n can be taken for a harmonic oscillator. In our case: 

n » 1, c n 2 ~ -7==exp - v __ . 6 
V27rn V 2n J 

Let us calculate now an average value of the oscillator coordinate x using the wave function 
(f^). Taking into account that only transitions n — > n ± 1 are important, one gets: 

< x(t) >^J2\ Cc n \ 2 cos((ujo + bn)t). (7) 

For small enough values of t we can replace the summation over n by integration. The result 
is Gaussian damping of amplitude of the oscillation according to: 

< x >~ exp(-nb 2 t 2 /2) = exp(-(t/r c ) 2 ), (8) 

where 

r" 1 = ( | 6 |) = u (Jn/2hu \v\)= uo^Ehuo/2 \ v |) (9) 
The periodicity of expression (|7|) gives immediately the revival period: 
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Note first of all that amplitude of the oscillations is proportional to yn. Thus our phemeno- 
logical theory gives the same amplitude dependence of collapse time as the theory of Ref. || . 
We however didn't use so far any information about microscopic properties of our system. 
It means that the effect is not related in our description to the fact of Bose condensation. 
Every system possessing excitations with long enough lifetime is suitable. The coefficient v 
can be calculated or measured and thus the theory contains no fitting parameters. 

According to Eq.(|9|) the collapse time r c decreases with increase of n. One must however 
take into account still that there is a restriction on the oscillation energy | v \ E « 1, or: 

1 << n « : (11) 

Tiujq | v | 



This inequality ensures applicability of the approximation of weak nonlinearity. It gives 



Uo^TiujQ I v | 

The effect can be observed if this r min is less that the oscillation damping due to dissipation 
which is always present. The theory we presented is a rigorous one at conditions (|TTJ) . 
To apply the theory to our problem of oscillations of Bose condensed gas in a harmonic 
magnetic trap one must calculate the nonlinearity coefficient v. According to Stringari ||, 
the low frequency oscillations with frequency co>o ~ ojh are well described by equations of 
hydrodynamics if the number of atoms N is large enough, (uh is a frequency of oscillation 
of an isolated atom in the trap.) To observe the effect under consideration here one must 
try to have as small a number of atoms as possible. We however will use the hydrodynamics 
approximation taking into account that it works well for the present experimental conditions. 

Calculation of the frequency shift (|1]) is straightforward but the procedure is a little bit 
cumbersome, because this shift appears at equations for third order terms with respect to the 
oscillation amplitudes. (Some results about amplitude dependence of the mode frequencies 
have been presented in Ref [ IT3 1 _) Having in view to estimate the order of magnitude of 



the effect we will restrict ourselves to the consideration of a spherical-symmetric breathing 
mode in an isotropic trap only. In this case the hydrodynamic velocity v has the form 
v(r) = u(t)r. Introducing a new variable w according to u — w/w one can reduce the 



equation for w to the form [1C]. \'\ I 



w + u 2 H w — LUtf/m 4 = . (13) 
Substitution w(t) = 1 + x(t) and expansion with respect to x gives: 

x + 5u 2 H x = 20u 2 H x 2 - 120u 2 H x 3 . (14) 

One has from (|T4"D: 

x ~ Acos((jj t),uj = V5u>h,5u>/ojo = (7/6)A 2 . (15) 



(In fact A is a fractional amplitude of the oscillations in the cloud radius.) 

The energy of oscillation can be calculated as twice the mean kinetic energy, E = 
J poV 2 dV, where po( r ) is the equilibrium density of condensate in a trap. Calculation gives 
E = (3/7)NfiA 2 , where /i is chemical potential of the gas, calculated in the Thomas-Fermi 
approximation. Thus the nonlinear shift coefficient v turns out to be: 
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This gives finally the collapse time for the mode under consideration: 



(16) 



18y/2 _j fiN ^0.79 UN 

Tc = ir 4 " 7m = w ° ~ V (17) 

It is reasonable to think that equation ( [PTD gives the correct order of magnitude for different 
collective modes, and holds even for a small number of atoms where the hydrodynamic 
approximation cannot be justified rigorously. 

Let us estimate r c assuming typical experimental conditions of the paper M: N = 4500, 



uj q /2ii = 264 Hz, A = 0.2. According to calculations of Ref. |T2| J^N/Tiloq = 104. Then 
r c = 246 ms. The authors of Ref. reported that the lifetime of the low m = mode is 110 
ms. The calculated collapse time is in reasonable order of magnitude agreement with the 
observed value, however both experiments and calculations ]ITJ demonstrate absence of 
the amplitude dependence of frequency for this mode. Thus damping of this mode cannot 
be explained by the present version of the theory. Experimental conditions of [0] are less 
favourable for observation of the collapse because of the larger number of the condensed 
atoms. 

One must take into account also that there is a different reason for nonlinear effects in our 
system because of trivial nonharmonic corrections to the magnetic trap confining potential 
which must be considered too. 

More accurate calculations are in progress now. I hope that these calculations will 
disclose a mode which will be most suitable for observation of this effect. I would like to 
stress here that the first step to confirm the presented scenario of the damping is to discover 
predicted amplitude dependence. 



It is reasonable to look for other suitable systems to apply the theory developed here. 
Liquid helium clusters and small piezoelectric samples are possible candidates. 

In conclusion a general phenomenological theory of mode collapse and revival for a system 
with weak nonlinearity is presented. The theory doesn't contain any fitting parameters. The 
collapse time is expressed through a classical nonlinearity coefficient of the system. This 
coefficient is calculated for the breathing mode of the Bose condensate in a isotropic trap. 
The collapse time turns out to be about 250 ms for typical experimental conditions. 

I would like to thank J. L. Birman for various enlightening discussions. 
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